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ABSTRACT 

Understanding  the  adsorption  of  water  vapor  onto  activated  carbons  is  important  for 
designing  processes  to  remove  dilute  contaminants  from  humid  gas  streams,  such  as 
providing  protection  against  chemical  warfare  agents  (CWAs),  or  against  toxic  industrial 
compounds  (TICs)  used  in  a  terrorist  chemical  attack.  Water  vapor  isotherms  for  Calgon 
BPL  granular  activated  carbon  (GAC),  military  C2A1  GAC,  electrospun  activated  carbon 
nanofibers  (ACnF),  Calgon  Zorflex™  activated  carbon  cloth,  and  Novaloid-based 
activated  carbon  fiber  cloth  (ACFC)  are  presented.  Of  particular  interest  are  the  ACFC 
isotherms,  which  exhibit  an  unusually  high  degree  of  hydrophobicity.  The  ACFC 
isotherms  also  show  a  correlation  between  water  vapor  adsorption  hysterisis  and  the  level 
of  activation.  Water  vapor  isothenn  models  from  the  literature  are  compared. 

INTRODUCTION 

Understanding  water  vapor  adsorption  onto  activated  carbons  is  important  for  adsorption 
processes  such  as  steam  regeneration  of  activated  carbons  or  the  removal  of  contaminants 
from  humid  gas  streams.  Water  can  compete  for  available  adsorption  sites  and  reduce  the 
capacity  for  the  carbon  to  capture  target  compounds.  Under  high  humidity,  water  can 
completely  fill  the  micropores,  nearly  eliminating  the  capacity  to  capture  any 
contaminants  more  weakly  adsorbed  than  water  (Huggahalli  and  Fair  1996).  The 
coadsorption  of  water  and  organics  onto  GAC  has  been  studied,  including  the  hysteresis 
effect  (Rudisill,  Hacskaylo  et  al.  1992). 

One  relevant  application  for  which  contaminants  are  removed  by  adsorption  from  humid 
gas  streams  is  providing  protection  against  a  chemical  attack  in  either  military  or 
homeland  defense  (terrorist)  situations.  Whether  purifying  air  with  a  respirator  for 
personal  protection,  or  purifying  air  with  a  filtration  system  for  collective  protection  of  a 
shelter  or  building,  outside  air  is  the  medium  to  be  filtered,  and  it  is  not  humidity- 
controlled.  Another  application  where  outside  air  may  be  filtered  is  in  large-scale 
painting  operations,  where  exhaust  volumes  can  range  to  the  thousands  of  cubic  meters 
per  second,  and  the  cost  of  humidity  control  is  prohibitive. 

Activated  carbon  fibers  (ACFs)  show  great  promise  as  a  new  adsorbent  and  catalyst 
support  in  competition  with  more  traditional  activated  carbons  (Mays  1999).  Carbon 


fibers  can  be  woven  into  cloth  and  then  activated  to  form  activated-carbon-fiber-cloth 
(ACFC).  ACFCs  are  generally  made  from  pitch,  rayon,  polyacrylonitrile  (PAN),  or 
phenolic  (Novoloid  (Hayes  Jr.  1981))  resin.  The  fibers  are  woven  into  a  fabric  and  then 
activated  in  steam  or  CO2  at  high  temperatures  (e.g.  800°C).  ACFCs  have  the  ability  to 
be  regenerated  by  direct  electrical  resistance  heating.  Resistance  heating  has  been  shown 
to  be  more  rapid  and  energy  efficient  method  of  regenerating  activated  carbon  than  with 
steam  or  inert  gas  (Petkovska  and  Mitrovic  1994).  Because  of  this,  ACFCs  have  been 
studied  recently  in  an  effort  to  develop  next-generation  energy-efficient  cyclic  adsorption 
systems  (Subrenat  and  Le  Cloirec  2006;  Vidal,  Ramirez  et  al.  2006;  Hashisho, 
Emamipour  et  al.  2007). 

The  reduced  organic  adsorption  in  humid  gas  streams  onto  ACFCs  has  been  studied  for 
both  equilibrium  (Cal,  Rood  et  al.  1996)  and  dynamic  (Sullivan,  Rood  et  al.  2001) 
adsorption.  Water  adsorption  isotherms  are  herein  presented  for  a  variety  of  ACFCs,  and 
Calgon  BPL  as  a  reference,  in  an  attempt  to  more  understand  the  adsorption  of  water  onto 
these  materials. 

EXPERIMENTAL  METHODOLOGY 

Carbon  selection  and  sample  preparation 

Carbons  examined  in  this  study  included  Calgon  BPL™  granular  activated  carbon 
(GAC),  military  C2A1  granular  activated  carbon  (GAC),  Activated  carbon  nano  fiber 
(ACnF),  Calgon  Zorflex™  activated  carbon  cloth,  and  Novoloid-based  activated  carbon 
fiber  cloth  (ACFC).  The  matrix  of  Novoloid-based  ACFCs  tested  consisted  of  four  levels 
of  activation  of  virgin,  oxidized,  and  ^-treated  fabrics. 

The  ACFC  used  was  American  Kynol  ACC-5092-10,  ACC-5092-15,  ACC-5092-20,  and 
ACC-5092-25.  The  last  two  digits  of  the  ACFC  part  number  are  related  to  the  level  of 
activation  and  are  analogous  to  the  nominal  surface  area  of  the  material,  i.e.  -10  is  ~ 

1000  m  /g,  -15  is  ~  1500  m  /g,  and  so  on.  The  Calgon  Zorflex™  material  tested  was 
ACC-FM10  and  ACC-FM 10-700.  The  ACnF  material  used  was  produced  by  eSpin 
Technologies  Inc.,  of  Chattanooga  TN.  The  ACnF  is  a  non-woven,  produced  by 
electrospinning  polyacetonitrile  (PAN),  and  then  activating  by  conventional  means. 

The  oxidized  samples  were  obtained  by  treating  virgin  ACFC  samples  (0.025  m  ,  4  to  7g) 
having  different  grade  of  activation  with  150  ml  of  1/1  (V/V)  concentrated  HNO3  and 
H2SO4  solution  for  5  days  (Dimotakis,  Cal  et  al.  1995;  Dimotakis,  Cal  et  al.  1995). 
Samples  were  then  washed  with  distilled  water  until  neutral  PH  conditions.  Samples  were 
then  drained  from  water  and  then  dried  in  an  oven  at  150°C  in  presence  of  200sccm  of  N2 
flow. 

The  H2-treated  ACFC  samples  were  heated  at  950  °C  in  pure  Hydrogen  environment 
(Menendez,  Phillips  et  al.  1996;  Li,  Quinlivan  et  al.  2002).  Each  sample  was  put  in  a 
quartz  boat  inside  a  5  cm  outer  diameter  quartz  tube  fitted  inside  a  temperature  controlled 
tube  furnace  (Lindberg  furnace  model  54232;  temperature  controller  unit  59344).  The 
samples  were  heated  in  N2  for  2  hours  at  200°C  to  remove  chemisorbed  oxygen  and  then 


the  gas  was  switched  to  Hydrogen  (Zero  gas  grade)  at  200  seem.  The  samples  remained 
at  950°C  for  3  hours,  then  heating  was  stopped  and  the  samples  were  left  to  cool  down. 
When  the  temperature  decreased  below  100  °C,  the  gas  was  switched  to  N2. 

Samples  were  labeled  according  to  the  precursor  followed  by  a  suffix  referring  to  the  type 
of  treatment,  for  instance  -15V  is  not  treated  (virgin)  ACC-5092- 15,  -15A  is  acid  treated 
(oxidized)  ACC-5092-15,  -15H  is  hydrogen  treated  ACC-5092-15. 

A?  Isotherms,  surface  area  and  pore  size  distribution 

N2  adsorption  isotherms  were  generated  by  volumetric  measurements  using  ASAP  2010 
instrument  (Micromeritics,  Inc.).  Samples  were  degassed  for  more  than  8  hours  at  150°C 
under  vacuum  until  a  pressure  of  4  -  6  pmHg.  Acid-treated  samples  have  to  be  degassed 
for  longer  periods  compared  to  the  virgin  or  the  hydrogen-treated  samples  to  reach  this 
degassing  pressure  limit.  The  adsorption  isotherms  were  obtained  at  a  relative  pressure  of 
10'6  to  0.99.  N2  adsorption  and  desorption  occurred  at  the  N2  boiling  point  (77K). 

The  standard  Brunauer,  Emmett,  and  Teller  (BET)  method  was  used  to  fit  the  adsorption 
data  over  a  relative  pressure  range  of  0.06  to  0.20.  Micropore  volume  and  pore  size 
distribution  were  detennined  using  the  3D  model  (Sun,  Chen  et  al.  1998).  The  model  is 
based  on  the  Dubinin-  Radushkevich  (DR)  equation  and  is  used  to  model  the  adsorption 
of  gases  on  microporous  adsorbents  with  slit  shaped  pores  between  4  A  and  100  A.  The 
average  micropore  width  was  detennined  based  on  pore  size  distribution  for  sizes  <  20  A 
while  the  total  pore  volume  was  detennined  based  on  the  volume  of  nitrogen  adsorbed  at 

-5 

saturation  and  using  Nitrogen  liquid  density  (0.808  cm  /g). 

X-ray  photoelectron  spectroscopy  ( XPS ) 

The  elemental  content  (C,  O,  N,  and  S)  at  the  surface  of  the  virgin  and  chemically  treated 
ACFC  samples  was  determined  using  XPS  technique.  The  analysis  was  conducted  using 
a  PHI  5400  (Perkin-Elmer,  Physical  Electronics  Inc.)  instrument.  A  spot  size  of  1  mm2 
was  targeted  with  Mg  Ka  radiation  or  Al  Ka  radiations  at  a  power  of  300  W  at  15  kV 
under  ultrahigh  vacuum  (10‘8  to  10'9  Torr). 

IhO  adsorption  isotherms 

H20  adsorption/desorption  isotherms  were  collected  gravimetrically  with  a  Gravimetric 
Sorption  Analyzer,  Model  GHP-F,  manufactured  by  VTI  Corporation  of  Hialeah,  FL 
(Figure  1).  The  system  uses  a  Magnetic  Suspension  Balance  by  Rubotherm  of  Bochum, 
Germany.  The  matched  mass  flow  controllers  in  the  system  are  rated  for  a  nominal  1 .0 
sLpm  flow.  Carbon  samples  used  for  the  isotherms  ranged  from  40  mg  to  150  mg. 

Nitrogen  was  used  as  the  carrier  gas  throughout  the  experiments.  Each  sample  (except  the 
C2A1)  was  dried  at  250°  C  under  nitrogen  prior  to  each  isothenn  until  a  static  weight  was 
achieved  (<0.03%  change  in  10  minutes)  or  for  a  maximum  time  of  one  hour.  The  C2A1 
carbon  was  dried  at  105°  C  in  order  to  minimize  impact  on  the  impregnants  in  the  carbon. 

The  GPH-F  contains  algorithms  for  generating  a  specified  relative  humidity  (RH)  based 
on  the  evaporator  temperature,  the  sample  temperature,  the  flow  rates  of  the  mass  flow 
controllers,  and  a  correlation  for  the  vapor  pressure  of  water.  In  order  to  verify  the 


accuracy  of  the  RH  generated  by  the  instrument,  an  Edgetech  Vigilant™  Chilled-Mirror 
Hygrometer,  with  a  model  S3  sensor,  was  connected  in  series  with  the  GPH-F.  The  set- 
point  and  measured  RH  values  were  within  an  average  1%  RH.  The  total  gas  flow 
through  the  system  was  monitored  with  a  mass  flow  meter,  Omega  FVL-1600  0-5  LPM. 


CARRIER  GAS  SAMPLE  GAS 


Figure  1.  Gravimetric  Sorption  Analyzer  model  GPH-F  schematic  and  photo. 


EQUILIBRIUM  ISOTHERM  MODELING 

A  large  collection  of  mathematical  fonnulas  have  been  proposed  by  previous 
investigators  to  model  isotherms  of  adsorption.  Here,  we  present  an  overview  of 
equations  that  have  been  proposed  to  model  adsorption  of  water  onto  carbon.  Many  of  the 
more  recently  developed  models  have  only  been  fit  to  a  few  sets  of  experimental  data. 
While  our  approach  to  these  formulas  is  genealogical,  we  will  lit  all  these  equations  to 
our  data  to  examine  which  equations  better  model  the  spectrum  of  materials. 

We  will  attempt  to  elucidate  the  evolution  from  one  formula  to  another  and  briefly  state 
similarities  of  form  to  other  equations.  To  make  these  similarities  more  clear,  we  have 
sometimes  stated  the  equations  in  a  different  form  from  those  given  by  the  original 
authors.  Quite  a  few  recent  equations  take  the  form  of  a  summation  of  two  terms;  these 
include  the  Stoeckli  equation,  all  equations  in  the  Do  &  Do  family,  and  the  extended 
CMMS  equation.  Oftentimes,  when  fit  to  a  set  of  experimental  data,  one  of  the  factors 
preceding  the  term  representing  how  much  of  the  adsorption  that  temi  is  responsible  for 


(represented  by  a  with  a  subscript)  will  become  zero  or  close  to  zero;  this  indicates  that 
the  term’s  contribution  is  insignificant,  and  overfitting  is  occuring.  Usually  if  this  factor 
is  set  to  zero,  the  equation  will  reduce  to  a  previously-developed  isothenn  formula. 

It  is  important  to  evaluate  these  models  not  only  on  a  quantitative  basis  -  that  of  measures 
of  error  and  correlation  -  but  also  a  qualitative  basis.  Mathematical  qualities  that  should 
be  satisfied  including  having  a  finite  (usually  nonzero)  Henry’s  Law  limit,  some  certain 
fixed  values  ( a  =0  at  pr  =  0,  a=a0  at pr=  1),  and  being  increasing  over  the  entire  range  of 
RH  values.  They  should  also  be  independent  of  units  used.  It  is  also  worth  mentioning 
that  some  formulas  may  fit  to  experimental  data  better  more  due  to  the  flexibility  of  the 
form  of  their  equation  than  the  theory  used  to  derive  them.  If  a  formula  is  fit  to  a  data  set 
from  a  carbon  with  very  different  properties  than  the  assumptions  made  in  deriving  the 
model,  it  may  produce  a  good  fit  regardless,  though  the  fitting  parameters  may  not  have 
physical  meaning.  This  is  increasingly  common  as  the  number  of  parameters  to  a  formula 
becomes  larger.  Since  many  proposed  formulas  are  modifications  to  previous  formulas, 
we  will  group  them  and  examine  them  as  such. 

In  the  following  table,  pr  is  the  relative  pressure  of  water  (relative  humidity),  a  is  the 
amount  adsorbed  (in  mg/g),  T  is  temperature  (in  Kelvin),  Rg  is  the  universal  gas  constant, 
E0  is  a  reference  energy,  and  N  i s  a  fixed  value  indicating  how  many  terms  to  sum  over. 
All  other  values,  unless  defined  below,  are  treated  as  fitting  parameters;  except  for  <f,-  in 
the  Levan  equation,  all  fitting  parameters  must  be  positive. 
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DUBININ-RADU SHKE VICH  TYPE  FORMULAS 

The  term  ‘micropore  filling’  was  first  coined  in  Dubinin’s  adsorption  work;  since 
micropore  filling  is  important  to  these  carbons,  it  is  reasonable  to  start  there.  The  original 
Dubinin-Radushkevich  equation  dates  to  1947.  The  DR  equation  was  found  to  perfonn 
well  on  most  activated  carbon,  but  not  deal  well  with  carbons  with  very  fine  micropores. 
((Do  1998),  p.  1 56)  The  Dubinin- Astakhov  equation  replaces  the  exponent  of  2  with  n,  a 
parameter  that  increases  with  the  degree  of  activation. 

The  Stoeckli  equation,  created  to  fit  type  IV  isotherms  on  oxidized  carbon,  is  a  sum  of 
two  DA  equations,  one  to  account  for  the  type  1  behavior  and  one  to  account  for  the  type 
V  behavior. 

The  Doong-Yang  equation  is  one  of  only  two  equations  we  have  found  that  accounts  for 
hysteresis  in  desorption.  To  model  the  desorption,  hG  is  set  to  the  relative  pressure  at  the 
beginning  of  the  hysteresis  loop.  Gamma  is  an  empirical  activity  coefficient,  and  is  found 
using  the  Dubinin-Serpinsky  equation  DS2,  which  will  be  discussed  in  the  next  section. 

All  equations  of  this  family  have  a  Henry’s  Law  limit  of  0  (as  long  as  n>l)  ((Do  1998), 
p.162),  which  many  previous  investigators  have  remarked  is  thermodynamically 


inconsistent.  However,  in  very  hydrophobic  carbons,  the  value  of  the  Henry’s  Law 
constant  is  very  small,  and  this  limit  can  be  a  reasonable  approximation. 

DUBININ-SERPINSKY  TYPE  FORMULAS 

The  first  four  equations  of  DS  type  vary  in  the  factor  that  is  used  to  account  for  the 
decrease  of  available  adsorption  sites.  DS2  was  proposed  by  Dubinin  (Dubinin  and 
Serpinsky  1981),  and  DS4  by  Barton  et  al.  as  a  modification  to  the  DS2  (Barton,  Evans  et 
al.  1992).  There  are  members  of  the  DS  family  we  have  chosen  not  to  treat.  DS1  fits 
poorly  to  data  at  middle-to-high  RH  and  is  not  treated  here.  DS3  claims  that  the  relative 
humidity  when  a  =  0  is  usually  more  than  1/3,  and  that  the  derivative  da/dp  there  is  large; 
this  is  clearly  aphysical.  The  Corrected  DS  equation  (CDS)  was  proposed  by  Gauden  and 
contains  an  extended  empirical  factor  to  improve  on  DS4  at  relative  pressures  less  than 
10%,  but  due  to  the  limitations  of  our  experimental  system,  we  do  not  have  adsorption 
data  at  relative  pressures  between  0%  and  10%.  The  CDS  can  also  be  shown  to  not  be 
independent  of  units,  because  of  logarithms  and  exponentials  in  its  formula.  For  these 
reasons  we  have  chosen  not  to  treat  these  equations. 

The  LeVan  equation  is  a  very  flexible  formula  that  is  essentially  a  truncated  power  series 
of  n/p  in  terms  of  n.  ft  is  included  here  because  the  derivation  follows  that  of  the  DS 
equations,  and  it  contains  the  DS2  and  DS3  equations  as  special  cases  (N=  2  and  N=  3, 
respectively),  ft  is  also  evident  that  if  N  was  allowed  to  tend  to  infinity,  it  would  be  able 
to  contain  the  DS4  and  the  CDS  equations,  as  well  as,  in  theory,  any  other  isothenn 
formula.  Usually  AM  is  enough  to  get  a  good  fit.  Aside  from  E,0,  which  is  the  Henry’s 
Law  constant,  the  physical  meaning  of  other  parameters  to  the  equation  is  not  clear.  Also, 
because  of  its  flexibility,  care  must  be  taken  to  make  sure  that  its  parameters  are  such  that 
the  equation  is  increasing  and  positive  throughout  its  domain. 

The  stated  fonn  of  these  equations,  which  give  the  relative  humidity  in  terms  of  the 
loading,  is  the  opposite  of  the  usually  desired  form  and  is  awkward  to  use  in  practice. 

DS2  can  be  stated  to  give  the  loading  in  terms  of  the  RH  by  solving  a  quadratic,  but  this 
form  is  somewhat  unwieldy;  since  for  purposes  of  fitting  the  two  equations  are  equally 
useful,  we  have  chosen  to  use  the  stated  form. 

DO  &  DO  TYPE  FORMULAS 

The  equation  proposed  by  D.D.  and  H.D.  Do  to  model  water  isotherms  has  received  a 
good  deal  of  attention  since  its  proposal  in  2000.  One  of  Do  and  Do’s  key  assumptions  is 
that  water  inside  micropores  forms  clusters  of  molecules.  The  original  Do  &  Do  (DD) 
formula  assumes  clusters  of  5.  The  CIMF  isotherm,  proposed  shortly  afterwards, 
generalized  Do’s  equation  to  clusters  of  arbitrary  size,  and  Cossarutto’s  modification  to 
the  DD  formula  went  one  step  further  and  simplified  the  portion  of  the  equation 
representing  adsorption  on  functional  groups.  The  first  term  in  these  three  models  is  a 
SIPS  equation,  and  the  second  term  of  Cossarutto’s  equation  is  a  Langmuir  isotherm. 


Lagorsse  et  al.  proposed  the  second  formula  that  accounts  for  desorption  (Lagorsse, 
Campo  et  al.  2005).  They  follow  a  similar  derivation  to  the  DD  equation,  but  assume 


clusters  of  7  molecules.  The  first  term  of  both  equations  of  the  Lagorsse  fonnula  is  also  a 
SIPS  isotherm. 

Furmaniak  et  al.  proposed  a  “Heterogeneous  Do-Do  Model”  that  accounts  for  different 
energies  of  adsorption  at  different  functional  sites  (Furmaniak,  Gauden  et  al.  2005). 
However,  without  prior  knowledge  of  the  nature  of  functional  sites  on  the  material,  one 
must  make  guesses  as  to  the  number  and  properties  of  the  sites  the  material  may  have. 
Since  we,  unfortunately,  do  not  have  that  knowledge,  we  have  chosen  to  not  treat  this 
model  in  our  paper. 

We  note  that  the  summands  in  these  equations  have  explicit  forms  that  simplify 
computation. 

OTHER  FORMULAS 

The  Qi-Hay-Rood  formula  is  another  type  5  model,  derived  to  be  an  improvement  on  the 
DS4  equation.  The  QHR  is  nonzero  at pr=  0,  but  its  value  there  is  very  close  to  zero  (Qi, 
Hay  et  al.  1998). 

Mahle  proposes  a  type  5  formula,  based  on  the  assumption  that  the  pore  distribution 
approximates  a  Cauchy  distribution,  which  is  similar  in  shape  to  a  normal  distribution 
(Mahle  2002).  It  is  known  that  the  pore  distribution  in  activated  carbons  is  often  bimodal, 
so  this  assumption  may  not  always  hold.  The  general  method  for  deriving  an  isothenn 
from  a  distribution  of  pore  sizes  is  discussed  in  chapter  3. 10  of  Adsorption 
Analysis: Equilibria  and  Kinetics  (Do  1998). 

Talu  and  Meunier  derived  a  type  V  isotherm  that  assumes  water  molecules  are  more 
attracted  to  each  other  than  to  the  carbon,  and  therefore  once  water  is  adsorbed  to  a 
functional  group  other  waters  will  cluster  onto  that  original  water  molecule  (Talu  1996). 

Rutherford  proposed  a  model  based  on  cooperative  multi-molecular  sorption  theory 
(CMMS)  that  explains  well  the  effects  of  interaction  with  functional  sites,  but  does  not 
account  for  microporosity  or  condensation,  and  does  not  work  well  at  relative  humidities 
above  80%  (Rutherford  2003).  The  CMMS  model  can  account  for  type  II  and  type  V 
isotherms.  Assuming  a  type  V  isotherm,  the  CMMS  model  reduces  to  an  equation  called 
the  Ising  equation.  Rutherford  also  proposed  an  Extended  CMMS  model  to  better 
examine  carbons  with  minimal  micro-  or  nanoporosity  (Rutherford  2006).  The  form  of 
the  extended  CMMS  equation  is  a  sum  of  an  Ising  equation  and  a  Langmuir  equation. 

MOLECULAR  SIMULATIONS  OF  ADSORPTION 

With  increasing  amounts  of  computational  power  available  to  the  researcher,  molecular 
simulations  of  adsorption  have  become  more  popular.  Earlier  work  is  summarized  by 
Brennan  (Brennan,  Bandosz  et  al.  2001).  More  recently,  Liu  &  Monson  using  the  Gibbs 
Ensemble  Monte  Carlo  (GEMC)  method  produced  isotherms  for  a  carbon  assumed  to 
consist  of  slit  pores  with  no  active  sites,  and  postulated  that  a  microporous,  hydrophobic 
carbon  would  resemble  their  isotherms  (Liu  and  Monson  2005).  We  will  later  show  that 


the  isotherms  for  our  hydrogen-treated  Kynol™  material  bear  certain  similarities  to  these 
idealized  molecular  simulations. 
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Figure  2.  Water  vapor  isotherms  with  the  SPC  model  for  slit  pores  of  0.65  nm  (Liu  and 
Monson  2005).  Reproduced  with  permission. 

In  a  more  recent  paper,  Liu  and  Monson  preformed  Monte  Carlo  simulations  on  a  carbon 
made  up  of  graphite  platelets  with  randomly  distributed  active  sites.  The  isotherms  from 
those  simulations  strongly  resemble  the  isotherms  of  Calgon  BPL  and  similar  carbons 
(Liu  and  Monson  2006). 


RESULTS  AND  DISCUSSION 


XPS  of  Virgin,  and  Chemically  Treated  ACFC 

Oxidizing  the  ACF  enriches  the  ACFC  with  oxygen-containing  functional  groups  such  as 
carboxylic  acids,  quinones,  and  phenolic  hydroxides.  Treating  the  ACFC  samples  with 
the  mixture  of  concentrated  HNO3  and  FI2SO4  increased  the  average  oxygen  content 
among  the  samples  from  4.5%  to  23.9%  while  the  average  carbon  content  decreased  from 
94.7%  to  74.6%.  The  average  nitrogen  content  increased  from  0.5%  to  1.2%  while  the 
average  sulfur  content  increased  from  0.21  to  0.25%.  Since  the  change  in  the  sulfur  and 
nitrogen  content  of  the  oxidized  samples  was  small  and  disproportionate  to  the  increase 
in  the  oxygen  content,  this  indicates  that  the  increase  in  oxygen  content  is  not  due  to 
sulfate  or  nitrate  groups  from  FI2SO4  or  HNO3.  On  the  other  hand,  heating  the  ACFC 
samples  at  950°C  in  UHP  H2  decreased  the  average  oxygen  content  among  the  samples 
from  4.5%  to  1.5%,  the  average  nitrogen  content  from  0.5%  to  0.4%  while  the  average 
carbon  content  increased  from  94.7%  to  98.1%.  For  acid  and  hydrogen  treatments,  the 
change  in  oxygen  content  didn’t  seem  to  depend  on  the  activation  level  of  the  virgin 
ACFC  samples. 

Table  1.  ACFC  composition  as  determined  by  XPS 


Sample  name 

Carbon 

(%) 

Oxygen 

(%) 

Nitrogen 

(%) 

Sulfur  (%) 

ACFC-10V 

94.1 

5.4 

0.3 

0.3 

ACFC-15V 

95.3 

3.6 

1.0 

0.1 

ACFC -20  V 

94.2 

4.8 

0.6 

0.3 

ACFC-25V 

95.4 

4.3 

0.2 

0.2 

ACFC-10A 

75.1 

23.8 

0.9 

0.2 

ACFC-15A 

74.4 

24.3 

1.0 

0.2 

ACFC -20  A 

72.6 

25.3 

1.8 

0.3 

ACFC-25A 

76.4 

22.3 

1.2 

Not  available 

ACFC- 1  OH 

98.4 

1.2 

0.4 

Not  available 

ACFC-15H 

98.7 

1.1 

0.2 

Not  available 

ACFC -2  OH 

97.8 

1.6 

0.6 

Not  available 

ACFC-25H 

97.5 

2.1 

0.4 

Not  available 

N2  Adsorption,  Surface  Area,  and  Pore  Size  Distribution  of  Virgin  and  Chemically 
Treated  ACFC 

The  FI-treated  ACFC  depicted  comparable  BET  surface  area  to  virgin  ACFC.  Going  from 
ACFC10  to  ACFC25,  the  BET  surface  area  increased  for  both  the  virgin  and  the  FI- 
treated  ACFC.  On  the  other  hand,  BET  area  was  comparable  for  all  the  oxidized  ACFC, 
and  was  lower  than  that  of  virgin  ACFC.  The  BET  area  for  the  oxidized  ACFC  was  not 
related  to  the  level  of  activation  of  the  precursor  virgin  ACFC. 


As  is  the  case  for  BET  surface  area,  the  micropore  and  total  pore  volumes  for  the 
hydrogen  treated  samples  were  comparable  but  slightly  lower  than  that  for  the  virgin 
samples,  while  the  acid  treated  samples  depicted  much  lower  micropore  and  total  pore 
volumes  compared  to  the  virgin  samples.  In  terms  of  ratio  of  micropore  to  total  pore 
volume,  the  acid- treated  samples  depicted  a  reduction  in  microporosity  for  all  ACFC 
activation  levels.  The  hydrogen-treated  samples  manifested  reduction  microporosity, 
particularly  for  ACFC  1  OH  and  ACFC25H. 


Table  2.  Characterization  of  ACFC  samples  using  N?  adsorption 


Sample 

BET 

Surface 

(m2/g) 

Micropore 

volume 

(cm3/g) 

Total  pore 
volume 
(cm3/g) 

Microporosity 

(%) 

Average 
micropore 
size  (nm) 

Average 
pore  size 
(nm) 

ACFC  10- V 

849 

0.40 

0.40 

100.0 

0.67 

0.67 

ACFC15-V 

1335 

0.62 

0.62 

100.0 

0.76 

0.76 

ACFC20-V 

1566 

0.74 

0.75 

99.7 

0.84 

0.86 

ACFC25-V 

1763 

0.83 

0.83 

99.9 

0.88 

0.88 

ACFC  10-0 

793 

0.33 

0.33 

99.3 

0.70 

0.73 

ACFC  15-0 

994 

0.45 

0.46 

97.9 

0.82 

0.94 

ACFC20-O 

668 

0.29 

0.30 

95.0 

0.82 

1.08 

ACFC25-0 

873 

0.37 

0.40 

93.6 

0.85 

1.11 

ACFC10-H 

926 

0.39 

0.39 

99.9 

0.71 

0.71 

ACFC15-H 

1314 

0.60 

0.60 

100.0 

0.80 

0.80 

ACFC20-H 

1509 

0.70 

0.70 

100.0 

0.86 

0.86 

ACFC25-H 

1739 

0.77 

0.77 

100.0 

0.88 

0.88 

For  all  ACFC  grades,  acid  treatment  resulted  in  more  change  in  mesopore  size 
distribution  compared  to  hydrogen  treatment.  Acid  treatment  increased  the  mesopore 
size.  Such  change  resulted  in  increasing  the  average  pore  size  of  the  acid-treated  ACFC. 
On  the  other  hand,  hydrogen  treatment  doesn’t  seem  to  affect  the  mesopore  size 
distribution.  With  respect  to  the  micropore  size  distribution,  the  impact  of  acid  treatment 
was  most  pronounced  for  ACFC20  and  ACFC25  while  ACFC  10  and  ACFC  15  depicted 
smaller  change  in  both  micropore  and  mesopore  size  distributions.  For  all  ACFC  grades, 
micropore  size  was  much  larger  than  the  mesopore  size  reflecting  the  high  microporosity 
of  ACFC, 

While  the  average  micropore  size  was  less  influenced  by  the  chemical  treatment,  the 
average  pore  size  for  pores  between  4  and  100  A  was  much  larger  for  acid-treated.  This 
confirms  the  above  suggestion  that  acid  treatment  resulted  micropores  fusing  together  to 
form  larger  pores.  On  the  other  hand,  both  virgin  and  hydrogen-treated  samples  depicted 
similar  average  micropore  and  total  pore  sizes  indicating  small  influence  of  hydrogen  and 
heat  treatment  on  the  pore  structure  of  ACFC,  The  decrease  in  the  micropore  volume  and 
the  absence  of  prominent  effect  on  the  micropore  size  is  consistent  with  previous  results 
reported  by  other  researchers  (Mangun,  Benak  et  al.  1999). 


To  summarize,  hydrogen  treatment  has  very  small  impact  on  micropore  and  total  pore 
volumes,  average  micropore  size,  BET  surface  area.  On  the  other  hand,  acid  treatment 


decreases  the  micropore  and  total  pore  volumes,  microporosity,  BET  surface  area. 
However  acid  treatment  increases  the  average  pore  size  and  doesn’t  have  significant 
effect  on  the  micropore  size.  Hence,  these  results  suggest  that  the  addition  of  oxygenated 
functional  groups  by  treatment  with  HNO3  and  H2SO4  mixture  resulted  in  decreasing  the 
volume  accessible  to  N2.  It  might  be  that  the  functional  groups  have  blocked  the  pores  of 
the  ACFC.  The  large  micropores  which  are  readily  accessible  to  nitrates  and  sulfate  ions 
might  be  affected  more  by  pore  blockage  than  the  small  micropores  that  might  be  less 
accessible  to  those  ions.  This  can  justify  the  reduction  in  the  micropore  and  total  pore 
volumes.  The  increase  in  the  mesopore  size  might  be  explained  by  the  action  of  the 
concentrated  acid  in  etching  the  carbon  structure  and/or  by  fusion  of  large  micropores 
and  small  mesopore  forming  larger  pores. 

Water  Adsorption  Isotherm  Results 

Isotherm  plots  for  the  various  carbons  are  presented  in  Figures  3-7. 
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Figure  3.  Water  vapor  isotherms  at  25°  C  on  Calgon  BPF  and  C2A1. 


The  BPL  isotherm  is  in  good  agreement  with  previously  published  results  (Huggahalli 
and  Fair  1996)  (Rudisill,  Flacskaylo  et  al.  1992).  The  C2A1  carbon  is  BPL  with 
impregnants  added  to  enhance  adsorption  of  chemical  warfare  agents.  The  impregnants 
add  primary  adsorption  sites  and  reduce  the  total  pore  volume.  Thus  the  onset  of  water 
adsorption  for  the  C2A1  is  earlier,  and  the  total  water  adsorption  is  reduced  compared  to 
the  BPL. 


The  ACnF  exhibits  minimal  hysteresis  and  is  the  most  hydrophilic  of  all  the  carbons 
tested.  Because  the  precursor  of  this  material  is  PAN,  the  material  likely  has  more 
primary  adsorption  sites  than  the  non-nitrogen  containing  carbons.  The  XPS  analysis  of 
this  material  is  pending. 


The  Calgon  Zorflex™  FM10  and  FM 10-700  use  rayon  as  a  precursor.  Their  curves  most 
closely  resemble  those  of  BPL,  whose  starting  material  also  contains  a  significant 
quantity  of  cellulose  structures. 


Figure  5.  Water  vapor  isotherms  at  25°  C  on  ACC-5092- 10,  -15,  -20,  and  -25 

The  Kynol™  ACFCs  exhibit  a  strong  type  5  isotherm  according  to  IUPAC  classification. 
This  is  characterized  by  low  adsorption  capacity  at  low  RH  followed  by  a  sharp  rise  and 
then  a  plateau  with  little  increase  in  adsorption  capacity  with  increasing  relative 
humidity.  The  degree  of  hysteresis  increases  with  the  level  of  activation.  The  shape  of 
the  ACFC  isotherms  most  closely  resemble  those  of  pitch-based  carbons  examined  by 
Kaneko,  et.  al  (Kimura,  Kanoh  et  al.  2004;  Ohba  and  Kaneko  2007).  This  is  expected,  as 
the  pitch  and  phenolic  resin  carbons  are  chemically  similar,  as  in  the  case  of  the  BPL  and 
Zorflex™.  The  Kynol™  is  a  slightly  more  homogenous  material,  and  so  the  highly 


activated  samples  exhibit  more  idealized  behavior,  specifically  the  separation  at  the 
bottom  of  the  hysteresis  loop  is  greater. 


Figure  6.  Water  vapor  isotherms  at  25°  C  on  acid-treated  ACC-5092-10,  -15,  -20,  and  - 

25 

The  oxidation  of  the  ACFCs  transformed  the  isotherms  from  type  5  IUPAC  classification 
to  type  4.  The  materials  were  drastically  modified  in  both  chemical  composition  and  pore 
structure  as  noted  in  Table  2,  and  the  isotherms  reflect  this.  With  >  20%  oxygen,  the 
primary  adsorption  sites  are  plentiful  and  water  adsorption  increases  with  increasing  RH 
throughout  the  isotherm.  Although  this  treatment  yields  increases  water  adsorption  in  high- 


humidity  situations,  which  may  be  undesirable,  it  has  utility  in  also  increasing  adsorption  of 
other  difficult-to-adsorb  compounds,  such  as  hydrides  (Mangun,  Benak  et  al.  1999). 


Figure  7.  Water  vapor  isotherms  at  25°  C  on  hydrogen-treated  ACC-5092-10,  -15,  -20, 

and  -25 


The  hydrogen-treated  ACFCs  have  increased  hydrophobicity.  The  -10  material  is  most 
radically  changed,  going  from  negligible  hysteresis  to  significant  hysteresis.  As  the  pore 
size  is  relatively  unchanged,  this  is  an  indication  that  the  surface  chemistry  is  key  to  the 
transformation.  The  -15  material  also  is  altered  significantly,  with  the  hysteresis  loop 
being  widened  at  the  base.  The  -20  and  -25  materials  show  a  less  dramatic  general 
widening  of  the  hysteresis  loop. 


When  Figures  2  and  7  are  compared,  it  is  apparent  that  the  hydrogen-treated  Kynol™ 
ACFC  is  more  similar  to  the  idealized  slit-shaped  pore  model  than  previously  studied 
carbons:  The  hysteresis  loop  is  wider  and  more  squared  at  the  bottom  of  the  loop.  The 
plateau  is  of  a  similar  shape  and  slope.  The  adsorption  loops  commence  adsorption 
around  an  RH  of  60%. 

The  Liu  and  Monson  model  predicts  that  if  there  are  no  primary  adsorption  sites, 
condensation  can  still  occur  in  the  pore  below  the  saturation  pressure  if  the  pore  width  is 
optimal  -  about  0.65  nm.  The  ACFCs  have  a  pore  size  distribution  that  is  not 
homogenous,  and  so  there  are  some  pores  at  this  optimal  size.  The  water  that  condenses 
in  the  optimal  pores  then  becomes  a  secondary  adsorption  site  for  further  water 
adsorption.  Thus,  significant  water  adsorption  will  generally  commence  at  a  RH  of 
-60%  in  all  supennicroporous  materials  regardless  of  the  surface  chemistry. 


EQUILIBRIUM  ISOTHERM  MODELING  RESULTS 

We  applied  a  least  squares  fitting  algorithm  to  the  equations  listed  above,  and  calculated 
the  r2  values  of  the  equations  against  the  experimental  data.  The  results  are  summarized 
in  table  3. 

Table  3.  The  regression  coefficient  r 2  for  each  pair  of  equation  and  experimental  data. 


DR 

DA 

Stoeckli 

DS2 

DS4 

LeVan 

DoDo 

CIMF 

Cossarutto 

ACnF 

95.53% 

99.18% 

99.61% 

98.23% 

99.85% 

99.39% 

99.81% 

99.81% 

99.93% 

C2A1 

Canister 

98.76% 

98.79% 

99.86% 

99.14% 

99.75% 

99.93% 

99.57% 

99.79% 

99.86% 

BPL 

99.49% 

99.83% 

99.78% 

99.36% 

93.59% 

99.31% 

98.25% 

99.95% 

99.88% 

Zorflex 

FM10 

99.65% 

99.70% 

99.90% 

98.86% 

99.37% 

99.42% 

99.83% 

99.88% 

99.88% 

Zorflex 

FM1 0-700 

99.14% 

99.70% 

99.62% 

99.32% 

99.31% 

99.31% 

99.78% 

99.79% 

99.79% 

ACFC 10V 

95.05% 

99.62% 

99.51% 

95.40% 

99.86% 

98.47% 

96.21% 

99.84% 

99.94% 

ACFC 15V 

94.75% 

99.60% 

99.49% 

97.91% 

99.46% 

98.41% 

92.69% 

99.79% 

99.90% 

ACFC20V 

94.54% 

99.77% 

99.71% 

98.28% 

99.58% 

98.81% 

91.56% 

99.90% 

99.93% 

ACFC25V 

93.13% 

99.95% 

99.94% 

97.78% 

89.25% 

98.42% 

88.79% 

99.93% 

99.93% 

ACFC10H 

92.82% 

99.89% 

99.97% 

97.47% 

98.79% 

98.02% 

88.37% 

99.94% 

99.94% 

ACFC15H 

94.26% 

99.94% 

99.93% 

98.68% 

93.37% 

98.77% 

88.39% 

99.95% 

99.96% 

ACFC20H 

98.26% 

99.77% 

99.70% 

99.54% 

99.52% 

99.52% 

93.49% 

99.81% 

99.82% 

ACFC25H 

98.15% 

99.73% 

99.69% 

98.10% 

89.93% 

98.13% 

92.09% 

99.80% 

99.81% 

ACFC15A 

99.36% 

99.33% 

99.84% 

99.27% 

98.98% 

99.22% 

99.87% 

99.85% 

99.89% 

Max 

99.65% 

99.95% 

99.97% 

99.54% 

99.86% 

99.93% 

99.87% 

99.95% 

99.96% 

Min 

92.82% 

98.79% 

99.49% 

95.40% 

89.25% 

98.02% 

88.37% 

99.79% 

99.79% 

Lagorsse 

CMMS 

full 

XCMMS 

CMMSV 

QHR 

Talu 

Mahle 

Max 

Min 

ACnF 

99.12% 

99.73% 

99.85% 

99.75% 

99.58% 

98.61% 

99.87% 

99.93% 

95.53% 

C2A1 

Canister 

98.71% 

99.45% 

99.87% 

99.49% 

99.46% 

98.92% 

99.76% 

99.93% 

98.71% 

BPL 

98.01% 

99.94% 

99.78% 

99.82% 

99.84% 

99.50% 

99.50% 

99.95% 

93.59% 

Zorflex 

FM10 

98.68% 

99.83% 

99.83% 

99.84% 

99.84% 

98.97% 

99.73% 

99.90% 

98.68% 

Zorflex 

FM1 0-700 

98.77% 

99.77% 

99.68% 

99.72% 

99.76% 

99.06% 

99.49% 

99.79% 

98.77% 

ACFC10V 

98.42% 

99.93% 

99.97% 

99.94% 

99.84% 

98.99% 

99.98% 

99.98% 

95.05% 

ACFC15V 

95.89% 

99.90% 

99.93% 

99.91% 

99.83% 

99.03% 

99.92% 

99.93% 

92.69% 

ACFC20V 

94.57% 

99.95% 

99.95% 

99.96% 

99.93% 

99.36% 

99.93% 

99.96% 

91.56% 

ACFC25V 

92.93% 

99.97% 

99.91% 

99.92% 

99.91% 

99.66% 

99.80% 

99.97% 

88.79% 

ACFC10H 

92.67% 

99.99% 

99.97% 

99.97% 

99.95% 

99.23% 

99.87% 

99.99% 

88.37% 

ACFC15H 

91.55% 

99.89% 

99.88% 

99.90% 

99.95% 

98.98% 

99.72% 

99.96% 

88.39% 

ACFC20H 

93.60% 

99.79% 

99.77% 

99.80% 

99.88% 

98.98% 

99.52% 

99.88% 

93.49% 

ACFC25H 

91.78% 

99.80% 

99.78% 

99.81% 

99.88% 

99.06% 

99.44% 

99.88% 

89.93% 

ACFC15A 

96.94% 

99.70% 

99.90% 

99.69% 

99.36% 

98.98% 

99.75% 

99.90% 

96.94% 

Max 

99.12% 

99.99% 

99.97% 

99.97% 

99.95% 

99.66% 

99.98% 

99.99% 

Min 

91.55% 

99.45% 

99.68% 

99.49% 

99.36% 

98.61% 

99.44% 

88.37% 

2 

Max  and  Min  are  the  maximum  and  minimum  r  value  along  each  row  or  column.  The 
Max/Max  and  Min/Min  cells  contain  the  maximum  and  minimum  r  value  for  the  entire 
table. 

Most  equations  seem  to  fit  well,  at  least  numerically,  to  the  hydrophilic  carbons.  It  can  be 
seen  from  Table  [number]  that  the  virgin  and  hydrogen- treated  Kynol  fabrics  have  overall 
worse  correlation  with  these  equations  than  the  other  adsorbent  materials.  There  are, 
however,  some  equations  that  do  well  on  both:  the  CMMS-type  equations  and  the 
Stoeckli,  CIMF,  QHR,  Talu-Meunier,  and  Mahle  equations  all  have  minimum  r  values 
of  greater  than  99%  for  this  data  set.  In  fact,  the  largest  r,  99.99%,  occurs  for  one  of 
these  hydrophobic  carbons  with  the  full  CMMS  equation. 

The  r  values  can  be  misleading,  though.  It  is  enlightening  to  look  at  the  graphs 
themselves  and  look  for  aberrations.  The  DR  and  Do  &  Do  equations  can  be  seen  to  have 
very  significant  overshoot  near  RH  =  100%.  This  is  occasionally  experienced  with  the 
CIMF,  and  full  CMMS  as  well.  The  DR  and  Do  &  Do  equations  also  fail  to  capture  the 
steepness  of  the  S-shaped  part  of  an  adsorption  curve.  The  LeVan  equation  sometimes 
returns  a  fit  that  is  not  monotonic,  such  as  in  ACFC10V.  The  Lagorsse  equation  does  not 
fit  very  well  on  hydrophobic  carbons  such  as  ACFC20FI,  and  for  hydrophilic  carbons 
such  as  ACFC15S  often  gives  a  desorption  branch  that  is  lower  than  the  adsorption 
branch,  which  is  an  aphysical  result.  Several  equations  show  strange  jumps  sometimes, 
such  as  Stoeckli  and  Do  &  Do  and  CIMF  on  eSpin  near  RFI  =  0%,  and  the  DS4  on  BPL. 

It  is  important  to  remember  that  these  correlations  are  for  these  specific  carbons, 
and  in  fact  most  of  the  equations  that  perfonn  poorly  here  usually  preformed  well  on  the 
carbons  they  were  originally  lit  to  in  their  original  papers. 
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Figure  8.  Water  vapor  adsorption  isotherms  at  25C  on  virgin  ACC-5092- 10,  with  plots  of 

fitted  equations. 


Figure  9.  Water  vapor  adsorption  isotherms  at  25C  on  acid-treated  ACC-5092-15,  with 

plots  of  fitted  equations. 


Figure  10.  Water  vapor  adsorption  and  desorption  isotherm  at  25C  on  acid-treated  ACC- 
5092- 15,  with  plots  of  fitted  equations. 
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Figure  11.  Water  vapor  adsorption  isotherms  at  25C  on  hydrogen-treated  ACC-5092-20, 

with  plots  of  fitted  equations. 


Figure  12.  Water  vapor  adsorption  and  desorption  isotherms  at  25C  on  hydrogen-treated 
ACC-5092-20,  with  plots  of  fitted  equations. 
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Figure  13.  Water  vapor  adsorption  isotherm  at  25C  on  Calgon  BPL,  with  plots  of  fitted 

equations. 
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Figure  14.  Water  vapor  desorption  isotherm  at  25C  on  Calgon  BPL,  with  plots  of  fitted 

equations. 


Figure  15.  Water  vapor  adsorption  and  desorption  isotherm  at  25C  on  Calgon  BPL,  with 

plots  of  fitted  equations. 


Figure  16.  Water  vapor  adsorption  and  desorption  isotherm  at  25C  on  C2A1,  with  plots 

of  fitted  equations. 
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Figure  17.  Water  vapor  adsorption  isotherm  at  25C  on  eSpin  ACnF,  with  plots  of  fitted 

equations. 


CONCLUSIONS 

There  is  a  variety  of  equilibrium  isotherm  models  in  the  literature.  These  models  were 
compared  over  a  diverse  group  of  activated  carbons.  Although  some  of  these  models  can 
fit  the  isotherms  with  an  r2  greater  than  99%,  none  of  these  models  are  predictive.  There 
are  no  isotherm  models  that  predict  hysteresis  with  any  reasonable  degree  of  accuracy. 

Observing  water  adsorption  isotherms  with  hysteresis  gives  qualitative  information  about 
the  chemistry  and  structure  of  the  carbons.  Carbons  that  are  similar  chemically  have 
similarly-shaped  isotherms,  e.g.  carbons  from  cellulosic  materials. 

Novaloid-based  ACFC  with  hydrogen  treatment  is  more  hydrophobic  and  has  a  more 
idealized  isotherm  than  any  previously  presented  activated  carbon.  This  material  would 
be  ideal  for  high-humidity  adsorption  applications. 

Although  no  real  material  can  be  modeled  as  simply  as  a  slit-shaped  pore,  the  ACFC  is  an 
ideal  material  for  the  subject  of  further  molecular  adsorption  simulations. 


NOMENCLATURE 


A 

Distribution  constant  for  Mahle  equation 

a 

Amount  adsorbed  at  a  given  pr 

ac 

Constant  for  DS4  equation 

af 

Capacity  of  functional  site  part  of  Do  &  Do  type  equations 

ai 

Capacity  of  type  I  term  of  Stoeckli  formula,  or  capacity  of  Ising  tenn  in 
Extended  CMMS  equation 

aL 

Capacity  of  Langmuir  term  of  Extended  CMMS  formula 

Capacity  constant  in  Talu  -  Meunier  equation 

Clo 

Pore  capacity 

Clp 

Concentration  of  primary  sites 

ay 

Capacity  of  type  V  term  in  Stoeckli  equation 

a n 

Capacity  of  microporous  part  of  Do  &  Do  type  equations 

B 

Distribution  constant  for  Mahle  equation 

c 

Henry’s  Law  coefficient  for  DS  equations 

E0 

Reference  energy 

H 

Proportionality  constant  in  Talu  -  Meunier  equation 

K 

Constant  for  Lagorsse  equation 

Ka 

Langumuir-type  constant  for  adsorption  part  of  Lagorsse  equation 

Kas 

Shape  constant  for  CMMS  equation 

Kd 

Langmuir-type  constant  for  desorption  part  of  Lagorsse  equation 

Kf 

Langmuir-type  constant  for  functional  site  part  of  Do  &  Do  type  equations 

K, 

Langmuir  type  constant  for  microporous  part  of  Do  &  Do  type  equations 

k 

Constant  for  QHR  equation 

k0 

Langmuir-type  constant  for  CMMS  equations 

ki 

Constant  for  CMMS  equations 

kL 

Langmuir  constant 

kt 

Constant  in  Talu  -  Meunier  equation 

N 

Number  of  terms  to  sum 

n 

Dubinin- Astakhov  exponent 

flj,  riy 

Stoeckli  exponents 

P50 

Relative  pressure  at  50%  of  a„  adsorbed 

Pr 

Relative  pressure  of  water  (relative  humidity) 

Rg 

Universal  gas  constant 

T 

Temperature 

P 

Energy  coefficient  of  an  adsorbate-adsorbent  system 

(PE0  =  E,  the  affinity  energy  of  the  system) 

Pi>  Pv 

Energy  coefficient  of  each  term  in  the  Stoeckli  formula 

K 

Fitting  parameter  for  Dubinin-Serpinsky  type  equations 

S 

Fitting  parameters  for  LeVan  equation 
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